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We evaluate the microwave admittance of a one-dimensional chain of fluxonium qubits coupled 
by shared inductors. Despite its simplicity, this system exhibits a rich phase diagram. A critical 
applied magnetic flux separates a homogeneous ground state from a phase with a ground state 
exhibiting inhomogeneous persistent currents. Depending on the parameters of the array, the phase 
transition may be a conventional continuous one, or of a commensurate-incommensurate nature. 
Furthermore, quantum fluctuations affect the transition and possibly lead to the presence of gapless 
“floating phases”. The signatures of the soft modes accompanying the transitions appear as a 
characteristic frequency dependence of the dissipative part of admittance. 
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I. INTRODUCTION 

Vortices of persistent current in superconductors have 
been viewed for a long time as a testing ground for 
various models of classical and quantum phase transi¬ 
tions. This is due to relatively strong interactions be¬ 
tween vortices coupled with a high degree of control 
over the vortex arrays. For example, inter-vortex inter¬ 
actions in the presence of a periodic external potential 
(created experimentally by modulation of the supercon¬ 
ducting film thickness 1 ) made vortices a convenient tar¬ 
get for investigation of commensurability transitions. 2,3 
Theoretical studies of the effects of randomness on vor¬ 
tex structure and vortex dynamics have led to the no¬ 
tion of collective pinning, 4,5 with importance stretching 
well beyond the physics of superconductivity. 6 Later, the 
discovery of high-temperature superconductors triggered 
studies of the vortex lattice melting transition and glassy 
behavior 1,8 in layered superconductors as well as struc¬ 
tural transitions of vortices interacting with columnar 
defects. 9 Because vortices in continuous superconductors 
have normal cores, vortex motion is dissipative 10 and 
much of the above-mentioned work used classical statisti¬ 
cal mechanics to address collective phenomena in vortex 
arrays. 11 

The effects of quantum fluctuations of vortices in 
continuous films become important only close to the 
supeconductor-insulator transition, which requires spe¬ 
cial tuning of the films’ normal-state resistance. 12 In ar¬ 
rays of Josephson junctions, however, vortices do not 
have cores, allowing for the study of quantum fluctua¬ 
tions. An array of small superconducting islands con¬ 
nected by Josephson junctions, where the island charg¬ 
ing energies were sufficiently large for quantum effects 
to be important, 13 was used in attempts to observe 
a quantum Kosterlitz-Thouless (KT) transition 14 and 
Mott transition 15 in a ID array. While these works 
advanced the nanofabrication techniques needed to pro¬ 
duce highly regular arrays of small Josephson junctions, 
the measurement results were ambiguous. The current- 


voltage ( I-V ) characteristics gave inconclusive evidence 
for a KT transition in single-line ID arrays. 14 ' 16,1 ' The 
observations of Ref. 15 were later interpreted to be consis¬ 
tent with a purely classical commensurability transition 
rather than the quantum Mott transition. 18 The suppres¬ 
sion of quantum effects in these experiments stemmed 
from the low inductance of the continuous superconduct¬ 
ing wires, which were necessary to make the Josephson 
junction arrays. 

We should note, also, that in the majority of experi¬ 
ments the evidence for the various classical and quantum 
transitions mentioned above was based on signatures in 
I-V characteristics. This method is limited to address¬ 
ing highly averaged quantities, and relies on substantial 
deviations of the investigated system from equilibrium 
(needed, for example, to overcome the static pinning of 
vortices). 

Recent developments in superconducting qubit tech¬ 
niques offer the possibility of overcoming the described 
limitations of previous experimental studies of many- 
body physics of vortices. Typical superconducting qubit 
experiments address the superconducting system using 
microwaves. 19 The system is only weakly perturbed by 
the microwave excitation; hence, this spectroscopic ap¬ 
proach allows one to probe the system close to equilib¬ 
rium. On the circuit element side, the development of 
the fluxonium qubit, 20 which combines a Josephson junc¬ 
tion with a superinductor 21 (i.e., an element exhibiting 
high inductance and low capacitance) opens avenues for 
studying quantum effects in superconducting arrays. 

In this work, we study theoretically a one-dimensional 
array consisting of superinductors and a chain of small 
Josephson junctions as shown in Fig. 1. Two parameters 
characterize such an array: the ratios of the Josephson 
(Ej) and inductive (El) energies (see Sec. II), combined 
into a characteristic length 

i = 2 \J Ej / E l , (1) 

and the ratio Ec/Ej of charging and Josephson energies, 
which controls the quantum fluctuations. Despite its sim- 
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plicity, this model allows for a variety of phase transitions 
of a classical or quantum nature (depending on the ratio 
Ec/Ej) as a function of the applied magnetic flux 4> e per 
plaquette of the array. 

Classically, for fixed £ > 1, there is a critical mag¬ 
netic field (f>* at which the system undergoes a transition 
from a homogeneous state with no persistent currents 
through the inductors to a state with static persistent 
currents. As a function of the characteristic length £, 
Eq. (1), we discuss two regimes, cf. Fig. 1(b): a “type-I” 
regime (1 < £ < y/2) featuring a second-order transition 
to a state of staggered persistent currents as shown in 
Fig. 2(a), and a “type-II” regime (£ 1), in which the 

transition takes place by the sequential (first-order) intro¬ 
duction of localized vortices of persistent currents in the 
lattice plaquettes, Fig. 2(b). These vortices correspond 
to kinks , meaning discontinuous jumps of height 27 t, in 
the node phase (or node flux) cj>j , see Fig. 1(a). Repul¬ 
sive interactions between kinks lead to a series of pinned 
commensurate phases with increasing (f> e . In the presence 
of quantum fluctuations, the number of kinks fluctuates 
and the initial transition at <fi* turns into a KT transi¬ 
tion. This is followed by commensurate-incommensurate 
transitions between classical pinned phases and quantum 
liquid phases of floating crystalline cells of a kink lattice. 

Each of the phases carries a “fingerprint” in the mi¬ 
crowave absorption spectrum, as the nature of low-energy 
excitations is sensitive to the types of phases and the 
transitions between them. The microwave spectra also 
carry information about the crossover from classical to 
quantum critical behavior in the vicinity of the transition. 
Yet another advantage of the spectroscopic approach is 
that it is a linear response to a weak perturbation. In this 
work, we will highlight the signatures of phase transitions 
that can be measured using microwave photons. 

The paper is organized as follows: In Sec. II, we formu¬ 
late the mathematical model for the circuit under con¬ 
sideration. We also introduce and discuss two methods, 
capacitative and inductive, for coupling the circuit to an 
external resonator, and give general forms for the radia¬ 
tion absorption rate in the two cases. 

In Sec. Ill, we describe the type-I regime, with £ com¬ 
parable to a lattice spacing. In this limit, the low energy 
excitations are gapped plasmon oscillations, where the 
node phases 4>j undergo small fluctuations. The plasmon 
excitations soften at the critical magnetic flux (f>*, lead¬ 
ing to a quantum Ising transition between a phase with 
the magnetic fluxes <pj = 0 (cf. Fig. 1) everywhere and a 
staggered phase <pj = (—1 ) J ^, with order parameter (j>. 

Section IV discusses the type-II {£ 1) regime. Here, 

the low-energy excitations are associated with the addi¬ 
tion or removal of kinks. As the magnetic field is in¬ 
creased from (f> e = 0, the cost of creating a kink is re¬ 
duced, leading to a visible peak in the absorption spec¬ 
trum below the plasmon continuum. At a critical field, 
this energy cost vanishes, leading to a proliferation of 
kinks and the formation of a series of gapped and gapless 
crystalline phases of kinks. Each of these has an observ- 



FIG. 1. (Color online) (a) Quantum circuit of coupled fluxo- 
nium qubits threaded by an external flux (f> e . (b) Phase dia¬ 
gram for the classical ground state, cf. Refs. 22,23 for the phase 
diagram of similar classical models. At inductances L such 
that 1 < l < y/2, the system is in the “type-I” regime that 
features a single Ising transition from a homogeneous phase 
into a phase of staggered persistent currents. For large L or 
£ 1 (“type-II” regime), the system, as a function of ex¬ 

ternal flux (j> e , develops subsequent phases of commensurate 
lattices of vortices or kinks. The rational numbers describe 
the kink density in the node flux configurations {cpj}. Darker 
shaded regions contain phases with higher denominators that 
complete the devil’s staircase. The line separating the homo¬ 
geneous and inhomogeneous phases is given by Eq. (9) in the 
type-I and by Eq. (30) in the type-II regime. 



FIG. 2. (Color online) (a) Staggered persistent currents char¬ 
acterize the high-field (cf> e > (j>V) ordered phase in the type-I 
regime, (b) A vortex of static persistent currents decaying 
over the length l, stable in the type-II regime. 

able signature in the excitation spectrum for adding or 
removing kinks. Finally, quantum effects such as broad¬ 
ening of peaks in the excitation spectrum and the ap¬ 
pearance of phases of incommensurate quantum liquids 
are discussed. 

A discussion of our analysis and its results is presented 
in Sec. V. 
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II. MODEL 
A. Lagrangian 

The circuit in Fig. 1(a) is described by the Lagrangian 

. 2 N—l 

£[</>, </>] = 2^7 (^' ~ 4-i) 2 - > ( 2 ) 

c j=o 

where = (2 e) 2 /C with C the Josephson junction ca¬ 
pacitance, and we formally set m = 0j- assuming pe¬ 
riodic boundary conditions. The potential in (2) takes 
the form 

JV-i F 

v M ~ Ej [ cos _ <k) _ !]} ( 3 ) 

i=o 

where 1?^ = $q/(47t 2 L), L is the inductance of the induc¬ 
tors and 4> e = 27t< 1>/d>o with $ denoting the magnetic flux 
per plaquette and $o = h/ 2e being the flux quantum. 
Furthermore, we define 9j = (j>j — 0j-i • the phase differ¬ 
ence across the Josephson junction between nodes j — 1 
and j, which in the following we denote as link j. 

The potential V[<f>] is not invariant under <pj —> cj)j-\-2n, 
as would be expected for a superconducting system. This 
is because we neglect the phase slip processes that al¬ 
low the inductors to relax to the true ground state of 
the system. Superinductors such as those used in a 
fluxonium qubit are engineered so that this is a valid 
approximation. 21 

At large C, the physics is dominated by the 
(meta)stable configurations {4>j} that minimize the static 
potential V[<p}. The magnetic flux <j> e and the characteris¬ 
tic length scale £, defined in Eq. (1), are the only param¬ 
eters in this “classical” regime. Throughout our theoreti¬ 
cal analysis, we assume the limit of large systems, £ -C TV. 
Note that, due to the inductive potential, a configura¬ 
tion with local flux 0 O is physically distinct from a state 
with flux (pj + 27m with integer n / 0. In the ground 
state, each variable cj)j will thus take values between — 7r 
and 7T. Furthermore, minimal configurations {4>j} neces¬ 
sarily satisfy ]Tb 0j = 0, which physically corresponds to 
current conservation. 


B. Response functions in absorption spectroscopy 

By coupling the circuit in Fig. 1(a) to microwaves, 
signatures of the various phases and phase transition 
become observable in absorption spectroscopy. Addi¬ 
tional interest in such experiments may arise because of 
well-defined peaks in the absorption spectrum below the 
plasma gap. We suggest two, in a certain way comple¬ 
mentary, schemes of coupling 19 the circuit to a microwave 
resonator: (a) inductive coupling and (b) capacitive cou¬ 
pling, see Fig. 3. 




FIG. 3. (Color online) Schematical antenna setups for absorp¬ 
tion spectroscopy: (a) capacitive coupling and (b) inductive 
coupling to a microwave resonator. 

1. Capacitive coupling 

In the situation of Fig. 3(a), we add an antenna capac- 
itively coupled to two neighboring nodes j a — 1 and j Q , 
which will couple the microwave mode to the Josephson 
junction at link j a . 

The Hamiltonian for the the microwave mode (with 
creation and annihilation operators b t and b) and its cou¬ 
pling to the circuit takes the form 

He = hujR^b + gc(b + b^)U ja . (4) 

Herein, uir is the frequency of the microwave mode and 
gc is the coupling constant, which is determined by 
the impedance of the microwave mode and the coupling 
capacitance. 19 Afj a = id/ddj a is the Cooper pair number 
operator for the antenna Josephson junction at link j a . 

Using Fermi’s golden rule, the absorption rate for 
driven radiation with frequencies uj incident on the an¬ 
tenna has the general form 

a c (u) =‘^^^\(m\Af ja \0}\ 2 5(u-LJ m ) (5) 

m 

with the sum being over the excited states | m) with en¬ 
ergy hujm relative to the ground state. 

2. Inductive coupling 

For the inductive scheme, Fig. 3(b), we add induc¬ 
tances La into the former “ground line” such that the 
system remains translationally symmetric. The anten¬ 
nas themselves may couple inductively to one or several 
of these inductances. Here we suppose it solely couples 
to one link, say j a - 

The inductances La should be chosen small such that 
the length 


£a = 2y/L/L A 


( 6 ) 


















4 


is not only large itself (in comparison with unity) but also 
l A £ (implying E La Ej), where £ is the character¬ 
istic length defined in Eq. (1). In this limit, as far as the 
ground state and excitation energies are concerned, the 
inductances La merely lead to a renormalized character¬ 
istic length £ e ff = £/y/\ + ( £/£a ) 2 - With good accuracy, 
we may thus set £ e ff — £ and ignore the effects of the 
inductances La in the analysis of the model (2). 

Using inductively coupled antennas, the microwave 
modes (described by b and tf) effectively couple to the 
Josephson phase differences 9j instead of Cooper pair 
numbers Mj, cf. Eq. (4). Currents induced by the an¬ 
tenna at link j a decay only over the (large) distance £a , 
Eq. (6). This leads to the following form of the absorp¬ 
tion rate: 

a L (uj) = IH0?JO)| 2 <5(w-w m ) (7) 


with 

^ _ 1 @k exp(ifcj a ) 

3a N k 1 + ^ sin2 (^/ 2 ) 


( 8 ) 


where k = 27m/AT with n = 0,1,..., N — 1 and 9k = 
9j exp(—i kj) is the Fourier transform of 9j. The last 
approximation is valid in the long-wavelength limit. The 
coupling constant gL in Eq. (7) is determined by the 
impedance of the microwave mode and the coupling in¬ 
ductances. 

In contrast to the capacitive coupling, which locally 
couples the microwave mode to one Josephson link, the 
inductive scheme has a much longer range £a , Eq. (6). 
We will specify the spectral absorption rates for the ca- 
pacitively coupled antenna, adai), and the inductively 
coupled antenna, «l(w), for the type-I regime in Sec. Ill 
and type-II regime in Sec. IV. In the following, we will 
use units in which H = 1. 


III. TYPE-I REGIME: MANIFESTATION OF A 
QUANTUM ISING TRANSITION 

For 1 < £ < y/2, as one increases the external flux <p e , 
the system undergoes a second-order phase transition 
from a low-flux homogeneous configuration <fij =0, cor¬ 
responding to zero persistent current through the induc¬ 
tors, into a staggered order <pj = <f>{— l)- 7 with an al¬ 
ternating persistent current through the inductors; see 
Fig. 2(a). This state is classically connected to the state 
with kink density p = 1/2 at larger £ in the type-II 
regime, cf. the phase diagram in Fig. 1(b). In this type-I 
regime, close to <j>*, the relevant low-energy excitations 
are plasmons, quantized small fluctuations in <pj. Mean- 
field theory, which assumes that these fluctuations are 
small, predicts that the transition happens at the critical 
flux <j)* given by 

cos<(>* = — 1/£ 2 for l<£<y/2. (9) 


Technically, it is more convenient to work with the 
field 9j = <pj — (pj-i ■ which represents the phase dif¬ 
ferences over the Josephson links. For <j> e > <f>*, the 
field 0, fluctuates around the ordered configuration 0, = 

20 (— 1 V- 


A. Mean-field theory 

For <j> e < (j>*, mean-field theory for the phase differ¬ 
ences 9j gives the plasmon spectrum: 


£ k = \VEcE^J^^:——+£‘ 2 cos(j) e , (10) 

2 V sin (fc/2) 

where, in the continuum approximation, the wave num¬ 
bers are k = 0,..., 27r. The lowest energy states occur 
for momentum k = n. For momenta k close to 7r, the 
spectrum is simplified to 


£fc “ 7 \/ 1 + l ^ (fc " 7r)2 ' (U) 

Herein, 

7 = \>\J E c E l ^/\ + £ 2 cos </) e 

~ y/E C Ej sin </)* |0 e -0*| 1/2 , </>e <</>:, (12) 

is the plasmon gap with the second line showing that 
mean-field theory predicts the gap closing as a square 
root as <j) e -A (j>*. At the transition, the softened low- 
energy plasmon modes have a linear dispersion relation 
£k = u\k — 7r| , where 



V EqEl 


(13) 


is the sound velocity. 

On the ordered side, for (f> e > </>*, plasmons describe 
the fluctuations of the order parameter <p. A similar 
mean-field theory calculation yields for momenta k with 
|fc — 7r| <C 1 a soft-mode dispersion relation of the same 
form as (11), with the gap larger by a factor of V2, 


7 = ~^=y/E C E L y/1 + £ 2 COS 4> e , 


> 


(14) 


but otherwise behaves as a function of <p e — 0* in the 
same way as the gap (12) on the disordered side. 


In absorption spectroscopy, microwave photons excite 
plasmons so that we expect manifestations of the critical 
point at <j) e = 4>* in the spectral absorption rate ac(«), 
Eq. (5), or ai(w), Eq. (7), for capacitive or inductive cou¬ 
pling to the microwave resonator. Specifically, we may 
write the Cooper pair number operator M :j and its con¬ 
jugate, the phase difference 9j = (pj — i n terms of 
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the plasmon modes, 


J\fj — 


M tak+al)e>k ’' 


d k l~Tk 

27r 


d k I Ec 


i=i y teV2s ( “‘ _ “‘ ) 


tw kj 


(15) 

(16) 


where a*, is the annihilation operator for a plasmon at 
wave number k. Inserting Eq. (15) into Eqs. (5), we find 
that at <j> e < (/>* the spectral absorption rate from the 
ground state is given by 

«cM = ^ ug(u>) (17) 

for the capacitively coupled antenna, where 

e(u) = J ^ S ( UJ ~ £fc ) ( 18 ) 

is the plasmon density of states. For inductive coupling, 
inserting Eq. (16) into Eq. (7), we find 


leading order in “discrete gradients” ftj — ftj-i so that 
16 + (fij ~ fij- 1 ) 2 ]- The Euclidean ac¬ 

tion then reads 


5 = 


Ej 


r/3/2 


dr 


-/3/2 


£{ 


(dr'&j 


E c Ej 4P 
+ -p.^3 - 2 COs[l?j - 0e(—1) J ] } 


{#j ~ #j- 1) 2 

(23) 


with /3 —> oo at zero temperature. Outside a Ginzburg 
region close to the critical field <)>*, one can treat the 
action (23) in the saddle-point approximation and the 
mean-field results from the preceding section become ac¬ 
curate. 

In order to determine the Ginzburg region, we ex¬ 
pand the cosine-potential to fourth order in ftj and 
then employ a continuum approximation, j —> x and 
ftj — ft j— i —t d x ft. Rescaling coordinates and fields so 
that they become dimensionless and so the prefactors of 
the quadratic terms are 1 / 2 , we find that close to <j>* the 
prefactor of the ft 4 term is small if 


a L (u>) 


ngjEc q(uj) 
i\ u 


(19) 


He 


> 


1 y/Eg/Ej 
12£ sin (j)* 


(24) 


for the softest plasmons with momentum k near 7 r, i.e. 
\k — 7 r| -C 1. Comparing Eqs. (17) and (19), we infer that 
the inductive-coupling scheme leads to stronger response 
at low energies w < Ec/£a and thus to higher-contrast 
results close to the transition at <f>*. 

The plasmon spectrum (11) implies the density of 
states 




1 cj&(uj — 7 ) 

TTU yjul _ q2 


( 20 ) 


where 0 denotes the Heaviside step function. In the fre¬ 
quency region W- 7 C 7 , the absorption rate features a 
square-root singularity, 


g(u) 


( 7 / 2) 1/2 6 ( u ; - 7 ) 

TTU y/UJ — 7 


( 21 ) 


At the critical field where the gap closes, 7 = 0, the 
van-Hove singularity disappears and Eq. (20) becomes 

Q{u) ^ — , (f> e = </* , (22) 

TTU 

as plasmons have become soft acoustic modes. 


B. Fluctuation regime 

Mean-field theory is valid as long as fluctuations are 
small. Close to the critical flux </>*, plasmon fluctuations 
become significant as the modes at k ~ 7 r soften. In¬ 
troducing ftj = (—1 yOj, we thus write an effective the¬ 
ory for the “slow” field ftj, keeping only the quadratic 


For external fluxes <p e satisfying Eq. (24), the mean-field 
results of the preceding sections are valid. 

Very close to <f>*, this Ginzburg criterion breaks down 
as quantum fluctuations become strong. The quantum 
critical behavior, which is due to the non-linearity of 
the Josephson current-phase relationship, leads to an 
excitation spectrum that is considerably different from 
mean-field theory and should correspond to a (1 + 1 )- 
dimensional quantum Ising chain . 24 As a result, e.g., the 
gap 7 in the plasmon spectrum, cf. Eq. (12), is expected 
to close at />* as 

'r~£t(E 3 jEo)*am{4>l) He~K\ > ( 25 ) 

i.e. with critical exponent v = 1 instead of v = 1/2 in 
the mean-field prediction (12). Furthermore, the t? 2 -term 
in the action (23) will become renormalized, effectively 
shifting the critical flux </>* to a higher value inside the 
Ginzburg region. At criticality, <p e = <f>*, the system 
is a liquid with spectrum = u(k — 7 r) and plasmon 
density of states of the form of Eq. (22), but in the fluc¬ 
tuation regime, the renormalized sound velocity u has 
to be considered an effective phenomenological param¬ 
eter. While parameters entering the observable quanti¬ 
ties ac(w) and ai(w), Eqs. (17) and (19), will be effec¬ 
tive ones, the qualitative threshold behavior should still 
be described in terms of Eqs. (21) and (22). 

Observability of quantum critical behavior requires the 
system size N to be larger than the correlation length at 
the boundaries of the quantum critical region as given by 
Eq. (24). This leads to the condition 

N^>r 1/2 (Ej/E c ) 1/4 ■ (26) 
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For typical parameters 25 for Ej and Eq (and £ ~ 1 in 
the type-I regime), the right-hand side of this estimate is 
of order unity. 

IV. TYPE-II REGIME: KINKS AND QUANTUM 
PHASE SLIPS 

Here we consider the limit £ P$> 1, corresponding to 
large inductances L. This regime is realizable with su¬ 
perinductors as demonstrated in fluxonium qubits. 20,21 
The most interesting physical effects are due to the pro¬ 
liferation of kinks in the node phases cj)j corresponding 
to current vortices as shown in Fig. 2(b). These vor¬ 
tices pick up currents over many plaquettes of the order 
of £ and therefore are stable already at small external 
fluxes <p e . Thus, phase transitions associated with kink 
proliferation preempt the instability driven by fluctua¬ 
tions of plasmons, which remain gapped for all external 
magnetic fluxes </> e . As a result, the ground state and 
excitations are fundamentally different from the type-I 
regime. 

We begin with the study of the phase diagram and the 
absorption spectrum in the “classical” limit. By “clas¬ 
sical” , we mean that the capacitances are large enough 
(and hence Ec is small enough) to make effects due to 
spontaneous quantum phase slips negligible but still al¬ 
low for induced phase slips by microwave absorption. 
Quantum effects due to finite capacitances alter the clas¬ 
sical picture and excitation spectrum and will be inves¬ 
tigated in Sec. IV C. 

A. Classical ground state 

The classical phase diagram in Fig. 1(b) is obtained 
by finding the configurations {(pj} of node fluxes that 
minimize the potential V[(j)\, Eq. (3), cf. Refs. 22 and 23. 
They are found from solving the set of equations given 

by 

1 dV 4 

° = Ejdfa = ~ [ sin ^' +1 “ ^ ~ ^ ( 27) 

- sin ((/)j - (j)j _i - 4> e )\ . 

Summation over j of Eq. (27) yields the con¬ 
straint 'pT J (pj = 0, which corresponds to zero net current 
to ground. 

At zero external flux cf> e , the ground state of the system 
is given by the homogeneous configuration <f>j = 0, which 
remains a local minimum of V[(j)\ for non-zero <f> e as long 
as cos (f>* > — l/£ 2 . In the limit of large £ 1, Eq. (27) 

admits non-trivial solutions already at small cf> e <C 1 
that contain kinks, 26 static local configurations of <f>j 
with a jump of the order of 2n across one link. In or¬ 
der to study such soliton solutions, we employ Villain’s 
approximation 27 , in which we expand sin(#j — </> e ) ~ 
0j — cf) e — 2nn with integer n such that 0j — cf> e — 2nn 
is small. 


For a single kink, with one jump of ~ 2n over the link 
between sites jo — 1 and jo and 9j = <p 3 — 4>j~\ < for 
all j ^ jo, Villain’s approach yields the configuration 

^ kink = -7rsgn(j + §) exp ( - 21j + \\/£) , (28) 

assuming that jo = 0. From Eq. (28), we see that £, 
Eq. (1), determines the scale of the width of a kink. In¬ 
serting typical fluxonium values, this width is £ ~ 6, 25 
which sets the smallest system size necessary to observe 
the physics under discussion. 

As the external flux (f> e is increased, we may expect 
the kink solution </> klnk , Eq. (28), to become more favor¬ 
able than the homogeneous configuration cj)j = 0 because 
of the Josephson junction’s preference for finite flux gra¬ 
dients at finite </> e , cf. Eq. (3). In fact, for the differ¬ 
ence A = V[<(> kmk ] — V[0] between the potential energies 
of the single kink and homogeneous configuration we find 

A = 2 nEj^l - <j> e ) , (29) 

where 

= j (30) 

is the critical external flux in the type-II regime, £ > 2, cf. 
Fig. 1(b). At flux (/>*, configurations {<j>j} with kinks be¬ 
come energetically favorable over the homogeneous con¬ 
figuration. The critical flux (j>* is analogous to the criti¬ 
cal magnetic field H c i in type-II superconductors, when 
vortices begin proliferating. 28 Its smallness in l/£ reflects 
that the vortex is able to pick up currents over the (large) 
length scale £, cf. Fig. 2(b). 

For <j) e > <t>*, the density of kinks in the ground state 
grows continuously as a function of <j) e . Each kink, when 
nucleating individually, brings an energy gain of A = 
2irEj((f>* — cf > e ). On the other hand, this individual gain 
has to be balanced with the interaction energy between 
two kinks. For the repulsive potential for two kinks at 
links i and j we find 

4? t 2 Ej ( 2\i — j \^ 

Jij = —j>— exp [ - j — ) , (31) 

which decays only at large distances on the scale of t, 
Eq. (1). 

The competition between kink-kink repulsion J t j, 
Eq. (31), and the single-kink energy A, Eq. (29), which 
plays the role of a chemical potential, completely de¬ 
termines the ground state as a function of the external 
flux <p e at £ 1. This interplay may be effectively de¬ 

scribed in terms of a classical spin chain model. Specifi¬ 
cally, we associate with each Josephson link j a pseudo¬ 
spin that distinguishes whether there is a kink across it 
(“spin-up”) or not (“spin-down”), see Fig. 4. The effec¬ 
tive spin Hamiltonian then has the form 

He l = A nj + — JijUiUj , 

j *7U 


( 32 ) 
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FIG. 4. (Color online) (a) Node flux (f>j configuration for a 
commensurate kink density p = 1/8 at £ — 3. (b) The same 
state represented in the effective spin ladder model for the 
coupled Josephson links. Links with upward orientied spin 
feature a kink in <f>j. 


where rij = (<rj + l)/2 and we introduce Pauli matri¬ 
ces crj, crj, and crj. 

The ground states of an Ising chain with infinite-range 
convex interaction potentials were systematically stud¬ 
ied as a function of the “chemical potential” A by Bak 
and Bruinsma 29 and Aubry 3 . In the language of our 
model, they showed that in the limit of very large sys¬ 
tems N —> oo, the kink density p = ( rij ) depends on cj) e in 
the form of a devil’s staircase: 30 This means that p only 
takes values in the rational numbers and rises monoton- 
ically and continuously with <f> e such that for each ratio¬ 
nal q/p < 1/2, there is a finite interval in </> e in which 
p{(j> e ) = q/p. For our effective model (32), this interval 
has the width 


A</> e (p) 


2irp 

£ 3 sinh 2 (p/^) 


(33) 


Figure 5(a) illustrates the devil’s staircase for fixed £, 
showing as a function of 4> e the kink density up to com- 
mensurability order p = 17. 

A ground state configuration {4>j} with kink den¬ 
sity p = q/p is periodic, 31 </>j+ p = </>y. For example the 
ground state configuration for p = 1/3 is • • • t|||||t|| 
• • •, with “primitive cell” til- For the non-unit fraction 
p = 2/5, a primitive cell contains a non-trivial basis: 

tltll- 

Figure 5(a) also illustrates that close to </>*, the depen¬ 
dence of the kink density on (f> e appears smooth so that 
we may try to effectively describe it using a continuum 
approximation. Physically such a description seems rea¬ 
sonable because close to 4 >*, kinks are sparse such that 
the large distances between neighboring kinks “wash out” 
the discrete structure of the underlying lattice. Thus, 
close to 4>l , the ground state density p of the classical 
model (32) is related to the external flux cf> e by 


Mp) 



1 

ip sinh 2 [l/(^p)]- 


(34) 



FIG. 5. (Color online) (a) Kink density p = q/p for the classi¬ 
cal ground state as a function of the external flux </ e (for i = 5 
fixed). The plot shows all densities with periodicity p < 17; 
higher-order commensurate phases are hidden behind the ver¬ 
tical dashed lines. The dashed curve is obtained in the con¬ 
tinuum approximation (34). (b) Classical excitation levels by 
adding (red, falling) or removing (blue, rising) a kink. For 
fixed 0 e , the lines mark the frequencies at which the spectral 
function exhibits a Apeak. The upper (lower) dashed curve 
shows the minimal energy for adding (removing) a kink in 
the continuum approximation, cf. Eqs. (44) and (52). For 
excitations close to </*, see Fig. 6. 


Solving this relation for p leads to the dashed curve in 
Fig 5(a). For (/> e — </>* </>*, Eq. (34) would predict the 

relation p ~ <j> e / 2-7T, which in particular implies p = 1/2 
at half flux quantum, cfi e = 7r, corresponding to the stag¬ 
gered order. However, discreteness effects are visibly su¬ 
perposed on the continuum model, mostly due to large in¬ 
tervals in (f> e that hold the same kink density, cf. Eq. (33). 

Very close to </>* such that ln[2</>*/(</> e — </>*)] 1, we 

explicitly find the dependence 



(35) 


which in particular shows that at cf> e —> </>* the kink den¬ 
sity grows with infinite slope. 

The transitions between the various commensurate 
phases of the model (32) are first order. 3 At zero tem¬ 
perature, quantum fluctuations due to the capacitative 
term in Eq. (2) facilitate equilibration and prevent hys¬ 
teresis effects if </> e is varied sufficiently slowly in an ex¬ 
periment. Realistic capacitive energies Ec ~ Ej lead 
to extended phases of incommensurate order, phases of 
“floating primitive cells”, around the transitions which 
we will discuss in Sec. IV C. 














































































B. Classical absorption spectrum 

In the regime of kink proliferation, the elementary ex¬ 
citations by photon absorption are the creation of an ad¬ 
ditional kink or the annihilation of a kink already in exis¬ 
tence in the ground state. The typical energy scale asso¬ 
ciated with such excitations is given by Ej/£ ~ ^EjEl- 
This scale could possibly already be below the plasmon 
gap ~ VE C Ej if El -C Eq- At the critical flux <f>*, 
Eq. (30), when the first kinks nucleate in the ground 
state, the spectrum for kink excitations extends to zero 
energy. Here, and in the vicinity of </>*, these are there¬ 
fore the relevant low-energy excitations. On either side 
of (j> e = </>*, the energies correspond to isolated low- 
frequency absorption peaks, see Fig. 6. 

In this section, we discuss the “classical” ab¬ 
sorption spectrum, i.e. we neglect spontaneous cre¬ 
ation/annihilation of kinks due to quantum phase slips 
and only allow for such processes in the context of pho¬ 
ton absorption. We discuss how quantum effects alter 
this picture in Sec. IV C. 

1. Capacitive vs. inductive coupling 

In a setup with a single capacitively- coupled antenna 
at link j a [Fig. 3(a)], the system may be excited by lo¬ 
cally introducing a kink at the link j a of the antenna. In 
a single absorption act associated with kink creation or 
annihilation at link j a , the system responds by boosting 
the phase difference 0j a over the j a -th Josephson junc¬ 
tion from 0 to 2ir or vice versa, whereas phase differences 
at other links acquire negligible changes < 1/1. For this 
reason, we may estimate the matrix element in the ab¬ 
sorption rate acH, Eq. (5), in the limit of a single fluxo- 
nium qubit 32 , where the energy difference s(j a ) between 
the excited and the ground state is determined by the 
external flux <fi e and the effective potential a kink at j a 
feels in the presence of kinks at other sites. 

The matrix element entering Eq. (5) is then given by 

\{e ja =2n\U ja \e ja =0)\^T/E c . (36) 

Herein, the parameter 

T = -^(E 3 jE c ) 1/4 exp ( - 8\/EjjEc) , (37) 

v 71 " 

with dimensions of energy, is the amplitude of a quantum 
phase slip. 33 In determining matrix elements and absorp¬ 
tion rates, we assume the typical limit Ej 3> Ec 3> 
El r and u> T for microwave frequencies to. 

As a result, we obtain for the absorption rate in the 
case of capacitive coupling 

a c (w) = 2 t rg% 8{u> - e(j a )) (38) 

with gc ~ gc • It features a single peak at fre¬ 
quency e(j a ), which corresponds to the energy cost for 


adding or removing a kink from the ground state at 
link j a The spectrum of energies e(j) is determined in 
Secs. IVB2 and IVB3, see also Fig. 5(b). 

The inductive coupling setup [Fig. 3(b)], as discussed 
in Sec. IIB 2, couples the single antenna to links over the 
long range of £a- For realistic situations, we may assume 
that l a is much larger than the length £, the relevant scale 
for the statics of kinks. This allows the insertion of kinks 
at any link over a long range. If the antenna is situated 
at j a , the absorption rate ai(w), Eq. (7), becomes 

Y^ e ~ 2b ~ 3a \ /eA $(u-£{j)) , 

(39) 

where <?l ~ <7l, as the matrix element in Eq. (7) is given 

by 

\{6j = 2Tr\dj\9 j = 0)| ~ r/e(j) . (40) 

For practical situations, because of the long range, it 
should therefore be sufficient to restrict ourselves to kink 
densities p = q/p with p -C l a- In this case, the quan¬ 
tum phase slip due to photon absorption may happen 
at any link within a ground state primitive cell, and the 
absorption rate is simply written as 

Qi ( u ) = 77 (w) [e+M + e-M] . ( 41 ) 

where g+(w) and g~(u) are the spectral functions associ¬ 
ated with the creation and annihilation, respectively, of 
a kink. We determine them in Secs. IV B 2 and IV B 3. 


2. Absorption-induced addition of kinks 

In this and the next section, we discuss the spectral 
function for excitations by adding or removing kinks. We 
focus mostly on the situation close to the critical flux cj>\. 
Figure 6 shows the photon absorption spectrum as a func¬ 
tion of the external magnetic flux (f> e . 

For external magnetic fluxes below the critical value 
<f>l, Eq. (30), the ground state is homogeneous, <f)j = 0, 
and photon absorption can only result in kink creation. 
The energy cost of a single kink is equal to A, Eq. (29), 
so that we find 

e+M = 5(to - 27tEj((/)* - 4> e )) for </> e < </>* • (42) 

The absorption frequency decreases linearly as a function 
of (f> e until it reaches 0 at the transition field </>*. 

For <f> e > (/>*, the ground state of the system carries a 
commensurate kink density p = q/p. Because of kink- 
kink repulsion, the energy costs for adding another kink 
at a particular site j is smallest if this link is put in the 
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e(u) 
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u/tP’s/EjEl 


0.5 



FIG. 6. (Color online) Classical absorption spectrum close 
to the critical flux </g. (a) Photon absorption spectrum as a 
function of external flux (j> e . Red solid lines correspond to ex¬ 
citations associated with the addition of a kink, blue dashed 
lines to removing a kink, (b.l) At </> e < </>* e , cf. Fig. 1(b), 
there is only one <S-peak in the spectral function q(oj), corre¬ 
sponding to introducing a kink into the homogeneous ground 
state, (b.2) For </ e > <$£, inhomogeneous broadening leads 
to a quasi-band of excitations related to kink creation, which 
for large energies turns into a series of (5-peaks. The lowest 
excitation is a single peak (for unit-fraction densities p = 1/p) 
corresponding to the annihilation of a kink. 


Thus, at the transition at (/>*, the gap (44) closes essen¬ 
tially as a square-root law. The absorption gap £+ in as 
a function of </> e is plotted as a dashed line in Fig. 5(b) 
and also as the lower band edge in Fig. 6(a). 

Adding a kink in the proximity of the center of the 
primitive cell is energetically least costly. Since (for p 
even) £+(_)) = e+(—j), each excitation level is doubly 
degenerate within one primitive cell. For small j—p/2 -C 
£, we expand Eq. (43) and find 


f 71 '—' £+ ■ 

^p\J ) — °mi 


8n 2 Ej(j~p/2) 2 

£ 3 sinh (p/£) 


(46) 


The accumulation of absorption levels close e)J) in turns 
into the usual van-Hove singularity in the continuum ap¬ 
proximation, which is valid for the low lying absorption 
levels close to </>*. Explicitly, for frequencies ui slightly 
above the energy gap, w — £+ in <C £^ in , the spectral func¬ 
tion for photon absorption associated with kink creation 
reads 


e+ M ^ 


s/t sinh(p/l) e(w- £+ in ) 
2tt{ P /£)V2E] 


(47) 


LO — £ ' 


for 4> e > (j>*. [We recall that by Eqs. (34) and (44) p is a 
function of </5 e .] Close to </>*, 


middle between two existing kinks, which are separated 
by a length of either the ceiling or floor integer of p -1 . 31 

Let us consider a magnetic flux slightly above (j )*. Here, 
kinks are sparse so that the discreteness of the lattice and 
peculiarities related to higher-order commensurate den¬ 
sities p = q/p with g / 1 are less important, at least 
for low-energy additional kinks. Since also, according 
to Eq. (33), most of the p!> e -space is filled by unit frac¬ 
tions p = 1/p, we focus at first on excitations by (low- 
energy) additional kinks to such ground states. With¬ 
out loss of generality, we may assume that the additional 
kink is introduced at link j = 0. Then assuming a ground 
state kink at j with j = 1,... ,p — 1 and for simplicity p 
even, the additional kink costs the energy 

+ («) 

Herein, the ground state density p = 1/p as a function 
of 4> e is found from inverting Eq. (34). For a given </> e , the 
lowest frequency at which a photon may create a kink is 
given by 


£ 


+ 

min 


£p{p/2) = A + 


47t 2 Ej 

lsmh(p/t) 


(44) 


which corresponds to adding a kink exactly in the center 
of the primitive cell of the kink lattice. For cf> e close 
above </>*, using Eq. (35), we find 


8 t t 2 Ej I 4> e - 4>* e 

Z y 2<f>% ln[2</>e/(<^e — </>*)] 


( 45 ) 


e+M 


V2 


Q (^ ~ £ min) 

\J (w/e+jn) - 1 


(48) 


with £ min S iven by Eq. (45). 

Equation (48) is valid for frequencies u> between e^ in 
and, roughly, 2e+ in . Integrating over this interval 
shows that the fraction of available kink-addition states 
with energies between e+ in and ~ 2e+ in is of order 
2/ln[2(/>e/(</> e — </>*)]. Close but not too close to </>*, this 
can already represent a significant fraction of all avail¬ 
able excited states associated with one additional kink. 
We also note that in the same limit close to <f>*, the level 
spacing |e+ (j) — £p(j —1)| remains smaller than the char¬ 
acteristic energy e+ in for the lowest (£/2)ln£ levels. In 
terms of energy, this corresponds to all states with ener¬ 
gies smaller than (H/2)e+ in . Both observations indicate 
that, as <p e —> <j)* the low-energy excitation spectrum as¬ 
sociated with the addition of a kink is suitably described 
by the continuum model (47) and (48) over a wide range 
as compared to the characteristic scale e+ in . 

At larger energies, however, the continuum approxi¬ 
mation underlying Eq. (47) eventually breaks down, even 
close to <j>*. Specifically, adding a kink right next to an 
existing kink determines a highest energy level 


£ 


+ 

max 


<(i) * 


4:TT 2 Ej 

T~ 


(49) 


The energy difference £+(l) — £+(2) between this high¬ 
est absorption level and the second highest, in which the 
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additional kink is placed on a next-nearest link, is of or¬ 
der Ej /£ 2 and to leading order independent of </> e — </>*, 
cf. Fig. 5(b). This invalidates any approach neglecting 
the discrete structure of the underlying lattice for finite 
l. The absorption spectrum at higher energies can then 
only be described in terms of a sequence of individual 
^-function peaks, 

0+M = -^25(uj-£+{j)) , (50) 

p zF 
3 

cf. Fig. 6(b.2). Note that the spectral weight of an in¬ 
dividual peak of Eq. (50) decreases linearly with 1/p, 
i.e. upon approaching the critical flux </>* from above, 
as 2/(£ln[2 <j)*/{(j) e — <(>*)]). In particular, the peaks in 
Eq. (50) are much less “bright” than the single peak on 
the homogeneous side, 4> e < <j>%, in Eq. (42). 

Let us finally discuss the situation of external fluxes 
away from the transition point, <j> e — <fi* > </*. For such 
external fluxes, the phase diagram for (j> e is dominated 
by phases with densities p = q/p where p is small. Ac¬ 
cording to Eq. (33), the phase with a ground state kink 
density of, e.g., p = 1/3 extends over an interval of length 
A</ e (3) ~ 27 t/ 3£, which is almost as large as the inter¬ 
val of the homogeneous phase at (f> e < (f>*. For finite £, 
this clearly invalidates the continuum approach entirely, 
and one can only work with the discrete formula for the 
absorption spectrum (50). As illustrated in Fig. 5(b), 
the absorption spectrum for the larger unit-fraction kink 
densities, which dominate the phase diagram on the </> e 
axis, contains fewer discrete 5-peaks, which on the other 
hand are individually brighter than the peaks close to </>*. 

Within an interval in <j> e of constant kink density, 
the excitation energies fall linearly as a function of <j> e , 
as a higher magnetic flux favors additional kinks, cf. 
Fig. 5(b). As a result, transitions between the commen¬ 
surate phases are detectable by jumps of the excitation 
frequencies seen in absorption spectroscopy. 


3. Absorption-induced annihilation of kinks 

For external magnetic fluxes above the critical flux, 
<f> e > <f>*, the ground state carries a finite density p of 
kinks whose annihilation constitutes another excitation 
in the context of photon absorption. In fact, excited 
states due to annihilation of a single kink are classically 
the energetically lowest excitations of the system. 

We again first consider a sparse unit-fraction den¬ 
sity p = 1/p at external flux slightly above </>*. In this 
situation, the energy cost to annihilate a kink is 


p £ exp(2 p/t) — 1 

Close to (f>*, using Eq. (35), we may write the dependence 


on <j> e explicitly, 

which is valid as long as \n[l<f)*/{(j) e — 1. Note that 

in contrast to the minimal energy to add a kink, Eq. (48), 
the energy for removing a kink depends on <f> e — <j>* in the 
leading order as a linear function. 

The spectral function then reduces to a single 5-peak, 

Q-(u>) = —8(ui - £~) for (j) e > <j>l . (53) 

P 

The spectral weight becomes small close to (j>* e , where p 
is large. This is in complete analogy with the fate of 
the individual 5-peaks in the spectral function £?+(w) for 
absorption-induced kink creation, cf. the discussion of 
Eq. (50). Observing this peak with a capacitively coupled 
antenna that affects the system locally thus seems much 
more difficult than with an inductively coupled antenna, 
which is able to excite remote links. 

For cj) e much larger than </>*, the ground state den¬ 
sity becomes larger and effects due to the discreteness of 
the lattice become important also for the spectral func¬ 
tion g_(u>). In particular, kink densities that are not unit 
fractions and feature more complicated primitive cells of 
the commensurate order may lead to a kink-annihilation 
spectrum with more than one peak. For instance, for 
the rational density p = 3/7, a primitive cell is given 
by 'M-'N-'N-I- Here, removing the second kink is energeti¬ 
cally less costly than removing the first or third, resulting 
in two 5-peaks in g_(u;) with the lower one being half as 
“bright” as the upper one, cf. Fig. 5(b) for different ex¬ 
amples of this phenomenon {p = 3/17 and p = 3/16). 

In conclusion of this section, the classical absorption 
spectrum shows clear signatures of the phase transitions 
between the various phases of commensurate kink config¬ 
urations and, in particular, the transition at </>*, at which 
the first kinks enter the system, cf. Fig. 6. The lowest ex¬ 
citations on the left side but essentially also on the right 
side of </>* are single 5-peaks. 

C. Quantum effects 

1. Quantum phase slips 

The capacitative term in the Lagrangian, Eq. (2), in¬ 
troduces quantum fluctuations into the system. Besides 
plasmons, which below the gap of order yJEcEj are 
frozen, these quantum fluctuations notably become man¬ 
ifest in the form of quantum phase slips, 3 which change 
phase differences across any of the Josephson junctions 
by 27 t. This corresponds to the spontanous creation or 
annihilation of kinks. For l 1, the necessary adjust¬ 
ments < £ _1 of fluxes 4>j neighboring the location of the 
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kink only make minor contributions to the action of the 
total quantum phase slip. 34 The phase slip amplitude T, 
which has already been introduced in Eq. (37), measures 
the coupling strength between the configuration with a 
kink at a given link and the one without. Including this 
coupling, the effective quantum Hamiltonian in the kink- 
dominated type-II regime £ 1 is 

h = h * + ^J2 <j * - ( 54 ) 

j 

where H c \ is the classical non-local Ising Hamiltonian, 
Eq. (32). The quantum phase slips thus enter in the 
form of an effective transverse field. As a result, in the 
fluxonium limit £ 3> 1, the circuit Fig. 1(a) constitutes 
a realization of a quantum Ising model with non-local 
interactions. 30 We remind the reader that the energy 
scales typically obey Ej 3> El 3> T. 


2. Low-frequency peaks 


The lowest excitations in the classical absorption spec¬ 
trum in the vicinity of </>* are single peaks, cf. Fig. 6 and 
Eq. (42). Quantum fluctuations broaden these peaks. 
Not too close to </>*, this effect is perturbative and the 
peaks should still be observable in the realistic quantum 
regime. 

For cf> e < (f >*, according to the classical analysis in 
Sec. IV B 2, the lowest excitation corresponds to intro¬ 
ducing a single kink into the homogeneous ground state. 
The excited state then has energy A, Eq. (29), and is 
IV-fold degenerate since the kink could have been excited 
on any of the N links in the system. As we include quan¬ 
tum effects due to the a x -term in the Hamiltonian (54), 
this kink, assumed at link i, effectively hops through the 
system: In second order in T, the kink either is first an¬ 
nihilated and then nucleates at a (different) link j or vice 
versa. In the first process, the effective hopping parame¬ 
ter is given by T 2 /A, in the second by — T 2 /(A + Vij) be¬ 
cause in the intermediate state the system holds two vir¬ 
tual kinks affected by the mutual repulsion Vij, Eq. (31). 
As a result, hopping is effective over distances 

£ = £\n(l/d), (55) 


where d, = A/(47r 2 Ej /£). The single-kink states thus 
reorganize into plane waves and the single excitation peak 
in Fig. 6(b) is broadened into a band. For |i — j | i, 
the two hopping processes interfere destructively. 

If A 2 IT 2 , we may use perturbation theory to find 
the effective quantum spectrum. To second order in T, 
the spectrum of excited states with a single kink is 


e n = A 


+ ^{*ln (1 + d - 1 )- 1 


(56) 


- tt£ 


sin[£g/2] 

svah(ir£q/2) 


+ {q O 2tt - q}') | 


with the wave number q = 0,..., 27t. In terms of <j> e , 
Eq. (56) remains valid as long as \<j> e — </>*\ £ 1 / 2 F/FJj. 

As (f> e —> <f>*, quantum fluctuations lead to drastically 
different low-energy physics. In the classical limit, T —> 0 
the spectrum reduces to a single peak at A, cf. Eq. (42). 

The bandwidth of the dispersive part in Eq. (56) is 
of order £T 2 /A. However, because the dispersive part is 
exponentially suppressed for momenta l/£ < q < 2i r — 
l/£, the band is mostly flat with most states accumulated 
at the lower band edge. For instance, a simple estimate 
using Eq. (56) shows that more than half of the states 
have an energy less than 


Ae 


^r 2 —7r 2 ^/4 

2A 


(57) 


away from e g=7r . The first term in braces in Eq. (56) 
represents a global renormalization oc T 2 of the classical 
single-kink excited level, due to interaction with virtually 
created and annihilated kinks. 


For cf) e > ()>*, the lowest energy single-kink excitation 
is the annihilation of an existing kink in the ground state 
(assumed unit-fraction for simplicity), cf. Fig. 6(b.2) and 
Eq. (53). At non-zero T, the induced hopping of kinks 
makes their positions subject to quantum uncertainty. 
As a result, the effective potential a single kink feels in 
the many-body background becomes fluctuating. The 
variance of this effective potential then determines the 
width of the broadened excitation peak for annihilation 
of a single kink. 

Perturbation theory in T yields a small broadening of 
the peak at e_, Eq. (52), with width 

Ae_ ~ £ 4 / 2 r . (58) 


Note the linear dependence on the quantum amplitude, 
in contrast to the broadening of the kink excitation spec¬ 
trum at (j> e < <t>*, cf. Eq. (57). The estimate (58) is valid 
for small T and for cf> e sufficiently far away from the crit¬ 
ical value (j)*, specifically as long as 


r< 


Ej (^ e - r e ) 1/2 

£ 3 / 2 ln 1//2 [2<^*/(0 e — <f>*)] ' 


(59) 


close to (j)*, and T -C Ej/£ 3 / 2 for cj) e n. If this in¬ 
equality is not satisfied, especially close to </>*, effects of 
the then strong effective hopping are more drastic and 
the single peak for kink annihilation fades away into the 
spectrum of a strongly-correlated quantum liquid, cf. the 
next section. Away from </)*, however, the peak should 
remain intact, owing to the exponential smallness of T, 
Eq. (37). 


3. Quantum phase diagram 

Quantum Ising models beyond nearest neighbors, i.e. 
Hamiltonians of the form of Eq. (54), have recently 35,36 
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FIG. 7. (Color online) Sketched quantum phase diagram, cf. 
Ref. 35, as a function of the external flux (p e for fixed l = 5. 
Quantum phase slips, in this illustration with amplitude F ~ 
0.1-Ej, melt all commensurate kink densities p = q/p with 
p > 6, leading to extended incommensurate Luttinger liquid 
(LL) phases, for which the spectral function has the form of 
Eq. (60). For general F, Eq. (61) determines the maximum p 
such that classical commensurate phases p = q/p survive the 
quantum fluctuations. 

been studied theoretically for power-law interactions, 
e.g., in the context of Rydberg atoms 35 and experimen¬ 
tally using trapped ions. 37 Applying the theoretical re¬ 
sults to our model (54), we predict the phase diagram as 
a function of external flux </> e for fixed i 1 and F to be 
the one shown in Fig. 7. 

Below a critical flux (/*, the system is, as in the clas¬ 
sical limit, in a homogeneous phase with gapped exci¬ 
tation spectrum, cf. Sec. IV C 2, but then undergoes a 
Kosterlitz-Thouless (KT) transition into a gapless phase 
of “floating” kinks. This phase is a Luttinger liquid and 
the excitation spectrum for adding or removing kinks is 
analogous to the spectrum of adding or removing spinless 
fermions in a one-dimensional system with strong repul¬ 
sion. Upon further increasing (f) ei the system undergoes 
a commensurate-incommensurate (CIC) transition into a 
“pinned” phase featuring a rational kink density p = q/p 
and a gapped excitation spectrum as in Fig. 5(b). Quan¬ 
tum effects broaden the “classical” peaks according to 
Eq. (56). After another CIC transition, the system be¬ 
comes again a Luttinger liquid. An (even) number of 
CIC transitions (depending on T) may follow before the 
system undergoes a second KT transition from a liquid 
into another gapped homogeneous phase, which in a final 
quantum Ising (QI) transition turns into the staggered 
order corresponding to kink density p = 1/2. 

Let us have a closer look into the phases that feature 
a finite kink density. For this purpose, let us consider 
a classical ground state with kink density p = q/p and 
estimate its stability to quantum fluctuations. Whereas 
the adjacent ground state for increasing (decreasing) </> e 
has only one extra (fewer) kink, practically all the kinks 
need to be rearranged in order to minimize the potential 
energy. 31 . We can interpret this state as the original state 
of period p with the addition of p defects, in analogy to 
the two domain walls introduced in an antiferromagnet 
by flipping one spin. 35,38 For the kink density p = 1/2, 
this analogy is perfect. 

The quantum a x term introduces effective hopping of 
the defects by virtual quasi-simultaneous adding and re¬ 
moving of kinks, similar to the single-kink hopping dis¬ 


cussed in Sec. IV C 2. The classical energies of the virtu¬ 
ally occupied intermediate states are of order Ej/t, cf. 
Fig. 5(b), giving a hopping amplitude t ~ T 2 £/Ej. The 
kinetic energy gain from such hopping ofp defects, which 
is ~ pt, then lowers the energy of the states with defects 
relative to the original period-p ground state and, as a 
result, takes away a strip of width ~ pt/2'KEj from both 
ends of the interval the period-p ground state has classi¬ 
cally occupied in <p e . In Fig. 7, e.g., the p = 1/5 phase 
had to cede almost one half of its classical interval in /> e , 
cf. Eq. (33), to the strips to its left and its right. 

In these strips, because of the mobility of the kinks, 
the commensurate configuration has melted and been re¬ 
placed by an incommensurate state of floating “defective” 
primitive cells. This state is gapless and physically de¬ 
scribed by a Luttinger liquid. Creation of a kink at a 
link j is here analogous to injecting a spinless electron at 
a site x into a (strongly-coupled) quantum wire. 39 There¬ 
fore, using the well-known Luttinger liquid results, we 
find for absorption at low energies the spectral function 

£±(w) oc oj^ k+k ~ 1 ~ 2) . (60) 

Within the incommensurate phase, the Luttinger param¬ 
eter K varies 36 as a function of <p e . At the CIC transition 
involving a commensurate kink density p = q/p, it takes 
the value K = 1/p 2 , cf. Refs. 35 and 39. If the width of 
the incommensurate strip exceeds the width of the inter¬ 
val of the classically commensurate state, i.e. if 

sinh(p/£) > , (61) 

the classical state with period p is entirely unstable and 
merges into an incommensurate phase of floating defects 
between two commensurate phases of shorter period p. 
The estimate (61) shows explicitly that classical ground 
states of kink densities with smaller denominator p are 
more immune to quantum fluctuations. In fact, the clas¬ 
sical phase p = 1/3, the “most immune” pinned phase, 
would be destroyed entirely only if T were so large that it 
would also destabilize the Luttinger liquid into a gapped 
homogeneous quantum paramagnet. 35 Based on typical 
fluxonium parameters 25 , we predict that classical phases 
of pinned densities p = q/p with p < 10 survive also in 
the presence the quantum fluctuations and systems of size 
larger than 10 would be needed to observe signatures of 
Luttinger liquid physics such as an absorption spectrum 
according to Eq. (60) . 

Upon decreasing of the external flux to a critical 
value (/>*, the lowest- </ e incommensurate phase destabi¬ 
lizes into a homogeneous phase in a KT transition 35 . 
We note that the value of the critical flux cj>*, because 
of quantum fluctuations, will be slightly larger than 
the classical value from Eq. (30). Close to this transi¬ 
tion, theory 35 predicts a critical K = K* = 1/8. As 
the external flux </> e approaches half a flux quantum, 
the system undergoes a similar KT transition from the 
rightmost (see Fig. 7) incommensurate phase (with the 
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same K*) into a homogeneous quantum paramagnet. 
Eventually, there is a quantum Ising transition 35,36 to 
the antiferromagnet-likc p = 2 phase in the vicinity of 

4>e = 7T. 

V. SUMMARY AND DISCUSSION 
A. Summary 

The model circuit, Fig. 1(a), we have discussed has a 
surprisingly rich equilibrium phase diagram, cf. Figs. 1(b) 
and 7, despite its relatively simple structure. This arises 
due to the combination of the nonlinear properties of 
Josephson junctions and the long-range interactions in¬ 
troduced by the coupling to a common ground. As 
a function of the circuit parameters l = 2^/Ej/El, 
Eq/Ej, and the external field (j> e , the model exhibits 
equilibrium phase transitions of the Kosterlitz-Thouless, 
commensurate-incommensurate, and Ising classes. We 
have shown that circuit QED realizations of this model 
enable access to low-energy excitations. In particular, 
they allow to identify the quantum phase transitions in 
linear response by the absorption of microwaves using 
a capacitively (C) or inductively (L) coupled antenna. 
Here we summarize our predictions for the absorption 
rate cxc/l(w) in the various characteristic regimes of pa¬ 
rameters, particularly l and <j> e . 

1. Type-1 regime: 1 < i < \[2 

In this limit, the elementary excitations of the system 
are plasmons. The absorption rates acM and ai(w) are 
given by Eq. (17) and (19), respectively, with the density 
of states of plasmon excitations given by Eq. (20). At the 
critical external flux <p e = <(>*, classically given by Eq. (9), 
the spectrum is gapless with a uniform density of states 
[Eq. (22)]. In the proximity of <f>*, where the Ginzburg 
criterion (24) is violated and quantum fluctuations are 
strong, the gap in the plasmon spectrum grows linearly as 
a function of the distance to </>*, [Eq. (25)], while outside 
the Ginzburg region, the classical square root dependence 
on the external flux [Eqs. (12) and (14)] sets in. 

2. Type-II regime: < > 1 

This regime, with Ej ^ El, is closer to the the pa¬ 
rameters realized in fluxonium qubits. 20 The ground state 
and excitation spectrum at large flux, cf> e > (f)* with <f>* 
given by Eq. (30), are significantly different from the 
small-flux regime, (f> e < <fr*, where the ground state is 
homogeneous. The elementary excitations are associated 
with the addition or removal of localized kinks (or vor¬ 
tices). 

The absorption spectrum is different for the capaci¬ 
tive [Eq. (38)] or inductive [Eq. (41)] coupling of the an¬ 


tenna, the former adding or removing kinks locally, the 
latter over an extended range. In the realistic quantum 
picture, the lowest excited states in the regime (j> e < </>* 
appear above a gap of order A oc (f>* e — (f> e and form a very 
flat quantum band [Eq. (56)]. The band slightly broad¬ 
ens the classical peak in the absorption rate [Eq. (42) and 
Fig. 6(b.l)]. As (j> e —> (/>*, the gap closes and, upon un¬ 
dergoing a Kosterlitz-Thouless (KT) transition at </>*, the 
system enters the phase of a “floating” crystal of kinks. 
The excitation spectrum here at <p e > <f>* is of Luttinger- 
liquid type and at low energies is given by Eq. (60). 

As (f> e is increased further, depending on the strength of 
quantum fluctuations (phase slips), the system undergoes 
a commensurate-incommensurate (CIC) transition, after 
which the ground state carries a pinned classical kink 
density. Here, the excitation spectrum is essentially the 
classical spectrum of Figs. 5(b) and 6(b.2), where quan¬ 
tum fluctuations slightly broaden the (5-peaks according 
to Eq. (58). The classical spectral function for kink ad¬ 
dition is given by Eq. (48), for (f> e not too far from </>*, 
and in general by Eq. (50). For kink annihilation, the 
spectral function is displayed in Eq. (53). For kink den¬ 
sities p = q/p with q = 1 or q = 2, which dominate the 
phase diagram, there is only a single peak in the absorp¬ 
tion spectrum related to kink annihilation, which is the 
lowest excitation energy of the system, cf. Eq. (51) and 
Figs. 5(b) and 6(b.2). 

As the external flux is further increased, classical and 
floating phases alternate as schematically indicated in 
Fig. 7. 


B. Discussion 

In our theoretical analysis, we built on the availabil¬ 
ity of a “superinductance” such as in fluxonium qubits 20 
when envisioning experimental realizations. As a result, 
the type-II regime of small El could be reached with¬ 
out considering the effects of additional parasitic capac¬ 
itances to ground that would be unavoidable in a real¬ 
ization using ordinary electromagnetic inductance. Such 
additional parasitic capacitances would suppress quan¬ 
tum fluctuations and thus enhance the classical behavior 
of the system. In the type-I regime, it decreases the ef¬ 
fective Ec for the low energy plasmon modes, narrowing 
the quantum critical region for the Ising transition. In 
the type-II regime, capacitance to ground decreases the 
quantum phase slip rate T, 34 which helps to stabilize the 
classical pinned phases. 

Throughout, we have used periodic boundary condi¬ 
tions for theoretical convenience. Realistically, it is eas¬ 
ier to create an array with open boundary conditions 
and so for small system sizes there will be edge effects. 
These will extend over a length ~ £, Eq. (1), as this is the 
scale for interactions in the system. Yet for large systems 
whose size significantly exceeds at least the length l, our 
results should be directly applicable. We have also as¬ 
sumed the temperature is zero. Both finite temperature 
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and the finite size of the system prevent the formation 
of long-range order and in principle mean that there will 
not be a true phase transition. However, the equilibrium 
state and spectrum will retain signatures of the infinite 
system behavior on short enough length and timescales. 

We note that there is an intermediate regime £ ~ 2, 
where the classical model has an incomplete staircase of 
first order transitions, cf. Fig. 1(b) and Refs. 22 and 23. 
In this parameter region, the type-II regime commensu¬ 
rate and incommensurate phases of finite kink density 
for <f> e > <f>* turn into the single Ising phase of the type- 
I regime. The classical theory 22 ’ 23 predicts a sequence 
of “superdegenerate” 22,23 and multicritical points as the 
transition at </>* changes from first order at large £ to sec¬ 
ond order in the type-I regime. The construction of an 
effective model for this region that also allows one to ana¬ 
lytically study the effect of quantum fluctuations appears 
difficult. However, the study of the ground states and ex¬ 
citations may be amenable to numerical techniques such 
as quantum Monte Carlo as the problem is bosonic. 

In the type-II regime at large £, the low-energy behav¬ 
ior is well described in terms of localized kinks. Using the 
locally-coupled capacitative antenna configuration, kinks 


can be selectively introduced or removed on individual 
sites. If multiple kinks are introduced, they will inter¬ 
act over the large length scale £ according the Hamilto¬ 
nian (32). This introduces the possibility of investigating 
the many-body non-equilibrium physics of the system. In 
the presence of dissipation, it may be possible to create 
a model system with a driven-dissipative steady-state of 
interacting kinks, an interesting addition to the set of 
non-equilbrium many body models that can be simulated 
with circuit QED systems. 40 


ACKNOWLEDGMENTS 

We would like to thank G. Rastelli, I. Pop, M. H. De- 
voret, J. Keeling, Yu Chen, C. Neill, S. Hacohen-Gourgy, 
and V. Ramasesh for interesting and helpful discussions. 
This work was supported by NSF DMR-1301798 and 
DMR-1206612, and ARO W911NF-14-1-0011. H. M. 
and R. T. B. acknowledge the Yale Prize Postdoctoral 
Fellowship. 


1 O. Daldini, P. Martinoli, J. L. Olsen, and G. Berner, 
Vortex-Line Pinning by Thickness Modulation of Supercon¬ 
ducting Films, Phys. Rev. Lett. 32, 218 (1974); P. Marti¬ 
noli, O. Daldini, C. Leemann, and B. Van den Brandt, 
Josephson Oscillation of a Moving Vortex Lattice, ibid. 
36,382 (1976); O. Daldini, C. Leemann, and P. Marti¬ 
noli, Coherent flux flow in superconducting films with peri¬ 
odically modulated thickness, Solid State Comm. 16, 509 
(1975); P. Martinoli, O. Daldini, C. Leemann, and E. 
Stocker, A.C. Quantum interference in superconducting 
films with periodically modulated thickness, ibid. 17, 205 
(1975). 

2 V. L. Pokrovsky and A. L. Talapov, Ground State, Spec¬ 
trum, and Phase Diagram of Two-Dimensional Incommen¬ 
surate Crystals, Phys. Rev. Lett. 42, 65 (1979) and The 
theory of two-dimensional incommensurate crystals, Sov. 
Phys. JETP 51, 134 (1980). 

3 S. Aubry, Exact models with a complete Devil’s staircase, 
J. Phys. C: Solid State Phys. 16, 2497 (1983). 

4 A. I. Larkin and Y. N. Ovchinnikov, Electrodynamics of 
inhomogeneous type-II superconductors, Sov. Phys. JETP 
38, 854 (1974). 

5 A. I. Larkin and Y. N. Ovchinnikov, Pinning in type II 
superconductors, J. Low. Temp. Phys. 34, 409 (1979). 

6 D. S. Fisher, Collective transport in random media: from 
superconductors to earthquakes, Phys. Rep. 301, 113 
(1998). 

7 L. I. Glazman and A. E. Koshelev, Thermal fluctuations 
and phase transitions in the vortex state of a layered su¬ 
perconductor, Phys. Rev. B 43, 2835 (1991). 

8 D. S. Fisher, M. P. A. Fisher, and D. A. Huse, Ther¬ 
mal fluctuations, quenched disorder, phase transitions, and 
transport in type-II superconductors, Phys. Rev. B 43, 130 
(1991). 


9 D. R. Nelson and V. M. Vinokur, Boson localization and 
pinning by correlated disorder in high-temperature super¬ 
conductors, Phys. Rev. Lett. 68, 2398 (1992) and Boson 
localization and correlated pinning of superconducting vor¬ 
tex arrays, Phys. Rev. B 48, 13060 (1993). 

10 J. Bardeen and M. J. Stephen, Theory of the Motion 
of Vortices in Superconductors, Phys. Rev. 140, A1197 
(1965). 

11 G. Blatter, M. V. Feigelman, V. B. Geshkenbein, A. I. 
Larkin, and V. M. Vinokur, Vortices in high-temperature 
superconductors, Rev. Mod. Phys. 66, 1125 (1994). 

12 A. M. Goldman and N. Markovic, Superconductor- 
Insulator Transitions in the Two-Dimensional Limit , Phys. 
Today 51, 39 (1998). 

13 R. Fazio and H. van der Zant, Quantum phase transitions 
and vortex dynamics in superconducting networks, Phys. 
Rep. 355, 235 (2001). 

14 E. Chow, P. Delsing, and D. B. Haviland, Length-Scale 
Dependence of the Superconductor-to-Insulator Quantum 
Phase Transition in One Dimension, Phys. Rev. Lett. 81, 
204 (1998). 

15 A. van Oudenaarden and J. E. Mooij, One-Dimensional 
Mott Insulator Formed by Quantum Vortices in Josephson 
Junction Arrays, Phys. Rev. Lett. 76, 4947 (1996); A. van 
Oudenaarden, B. van Leeuwen, M. P. M. Robbens, and J. 
E. Mooij, One-dimensional Mott localization of quantum 
vortices in Josephson-junction arrays, Phys. Rev. B 57, 
11684 (1998). 

16 D. B. Haviland, K. Andersson, and P. gren, Superconduct¬ 
ing and Insulating Behavior in One-Dimensional Joseph¬ 
son Junction Arrays, J. Low Temp. Phys. 118, 733749, 
(2000); K. Andersson and D. B. Haviland, Escape from a 
zero-current state in a one-dimensional array of Josephson 
junctions, Phys. Rev. B 67, 092507 (2003). 



15 


17 Y. Takahide, H. Miyazaki, and Y. Ootuka, Superconductor- 
insulator crossover in Josephson junction arrays due to 
reduction from two to one dimension, Phys. Rev. B, 73, 
224503 (2006). 

18 C. Bruder, L. I. Glazman, A. I. Larkin, J. E. Mooij, and A. 
van Oudenaarden, Phase transition in a chain of quantum 
vortices, Phys. Rev. B 59, 1383 (1999). 

19 A. Blais, R. S. Huang, A. Wallraff, S. M. Girvin, and R. 
J. Schoelkopf, Cavity quantum electrodynamics for super¬ 
conducting electrical circuits: An architecture for quantum 
computation, Phys. Rev. A 69, 062320 (2004). 

20 V. E. Manucharyan, J. Koch, L. I. Glazman, M. H. 
Devoret, Fluxonium: Single Cooper-Pair Circuit Free of 
Charge Offsets, Science 326, 113 (2009). 

21 N. A. Masluk, I. M. Pop, A. Kamal, Z. K. Minev, and 
M. H. Devoret, Microwave Characterization of Josephson 
Junction Arrays: Implementing a Low Loss Superinduc¬ 
tance, Phys. Rev. Lett. 109, 137002 (2012); M. T. Bell, 
I. A. Sadovskyy, L. B. Ioffe, A. Yu. Kitaev, and M. E. 
Gershenson Quantum Superinductor with Tunable Nonlin¬ 
earity, ibid. 109, 137003 (2012). 

22 C. S. O. Yokoi, L.-H. Tang, and W. Chou, Ground state 
of the one-dimensional chiral XY model in a field, Phys. 
Rev. B 37, 2173 (1988). 

23 M. Marchand, K. Hood, and A. Caille, Nonconvex in¬ 
teractions and the occurrence of modulated phases, Phys. 
Rev. Lett. 58, 1660 (1987) and Nonconvex interactions: A 
mechanism for the occurrence of modulated order in con¬ 
densed matter, Phys. Rev. B 37, 1898 (1988). 

24 See, e. g., S. Sachdev, Quantum Phase Transitions (Cam¬ 
bridge University Press, Cambridge, 2011). 

25 For estimates, we use the values Ej ss 10 GHz, El « 
1 GHz, and E c /8 = e 2 /(2C) « 3 GHz. 

26 M. Kardar, Josephson-junction ladders and quantum fluc¬ 
tuations Phys. Rev. B 33, 3125 (1986). 

27 J. Villain, Theory of one- and two-dimensional magnets 
with an easy magnetization plane, II. The planar, classical, 
two-dimensional magnet, J. Phys. 36, 581 (1975). 

28 P.-G. de Gennes, Superconductivity of Metals and Alloys 
(Benjamin, New York, 1966) 

29 P. Bak and R. Bruinsma, One-Dimensional Ising Model 
and the Complete Devil’s Staircase, Phys. Rev. Lett. 49, 


249 (1982). 

In order to avoid any confusion, we point out that although 
our interaction Jij, Eq. (31), falls exponentially and so is 
certainly short-ranged in the sense of scaling, it is never¬ 
theless non-zero at an arbitrarily large distance \i — j\. This 
implies the theoretical completeness of the devil’s staircase 
in our model, even though the “stairs” for kink densities 
p — q/p with large p occupy exponentially small intervals 
in (f> e , cf. Eq. (33). 

J. Hubbard, Generalized Wigner lattices in one dimen¬ 
sion and some applications to tetracyanoquinodimethane 
(TCNQ) salts, Phys. Rev. B 17, 494 (1978). 

G. Catelani, R. J. Schoelkopf, M. H. Devoret, and L. 
I. Glazman, Relaxation and frequency shifts induced by 
quasiparticles in superconducting qubits, Phys. Rev. B 84, 
064517 (2011). 

K. A. Matveev, A. I. Larkin, and L. I. Glazman, Persistent 
Current in Superconducting Nanorings, Phys. Rev. Lett. 
89, 096802 (2002). 

G. Rastelli, I. M. Pop, and F. W. J. Hckking, Quantum 
phase slips in Josephson junction rings, Phys. Rev. B 87, 
174513 (2013). 

E. Sela, M. Punk, and M. Garst, Dislocation-mediated 
melting of one-dimensional Rydberg crystals, Phys. Rev. 
B 84, 085434 (2011). 

P. Fendley, K. Sengupta, and S. Sachdev, Compet¬ 
ing density-wave orders in a one-dimensional hard-boson 
model, Phys. Rev. B 69, 075106 (2004). 

P. Richerme, C. Senko, S. Korenblit, J. Smith, A. Lee, R. 
Islam, W. C. Campbell, and C. Monroe, Quantum Catal¬ 
ysis of Magnetic Phase Transitions in a Quantum Simula¬ 
tor, Phys. Rev. Lett. Ill, 100506 (2013). 

A. A. Odintsov, One-dimensional Josephson arrays as su¬ 
perlattices for single Cooper pairs, Phys. Rev. B 54, 1228 
(1996). 

T. Giamarchi, Quantum Physics in One Dimension 
(Clarendon Press, Oxford, 2004). 

A. A. Houck, H. E. Tiireci, and J. Koch, On-chip quantum 
simulation with superconducting circuits, Nat. Phys. 8, 292 
( 2012 ). 


30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 


